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We show that magnètic susceptibility can reveal spin entanglement between individual con- 
stituents of a sòlid, while magnetization describes their local properties. We then show that magne- 
tization and its variance (equivalent to magnètic susceptibility for a wide class of systems) satisfy 
complementary relation in the quantum-mechanical sense. It describes sharing of (quantum) infor- 
mation in the sòlid between spin entanglement and local properties of its individual constituents. 
Magnètic susceptibility is shown to be a macroscopic (thermodynamical) spin entanglement witness 
that can be applied without complete knowledge of the specific model (Hamiltonian) of the sòlid. 
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Thermodynamical properties, such as heat capacity, magnetization or magnètic susceptibility, are normally ascribed 
to macroscopic objects with the number of individual constituent of the order of 10 23 . In contrast, genuine quantum 
features like quantum superposition or entanglement are generally not sccn bcyond molecular scalcs. As mass, size, 
complexity and/or temperature of systems increase the observability of their quantum effects is gradually limited by 
decoherence - an interaction of the system with its environment - that turns them into classical phcnomena. This raises 
several qüestions: under which conditions can quantum features of individual constituents of a sòlid have an effect on 
its global properties? Can one detect existence of quantum entanglement in a sòlid by observing its thermodynamical 
properties only? Can one consider macroscopic properties as quantum-mechanical observables in the sense that they 
obey complementary relations like position and momentum? 

The complementarity principle is the assertion that there exist observables which are mutually exclusive in the 
sense that they cannot be precisely defined simultaneously. One of them, for example, the z component of the 
spin- ^ (o~ z ), might be well defined at the expense of maximum uncertainty about the other orthogonal directions 
(a x and a y , <r's are respective Pauli matrices). One can speak about sharing of (quantum) information between 
mutually complementary observables Q. In the case of a qubit this can quantitatively be described by the relation 
(dj;) 2 + (<x y ) 2 + (er z ) 2 < 1, where the average is taken over an arbitrary state. When extended to composite systems 
the principle of complementarity asserts the mutual exclusiveness between entanglement and local properties of 
individual constituents of the composite system. In the case of two qubits this can be described by the relation 
Tli=x y z( a i) 2 + ( a i) 2 + ( a l a i) 2 — where the upper indices indicate qubits. The maximal value of 3 can be achieved, 
e.g., with product states (e.g. (al) = (<r 2 ) = (o^er 2 ) = 1; others are zero) for which local properties of the qubits are 
well-defined, but there is no entanglement. Alternatively, their joint properties could be well-defined at the expense 
of a complete mdefmiteness of the local properties (e.g. for a singlet state (o^er 2 ) = (o^cr 2 ) = (c^cr 2 ) = —1; others are 
zero) . 

Recently, a complementarity relation was proposed 0] between two macroscopic quantities, magnetization and 
magnètic susceptibility along one spatial direction. However, because entanglement between spin-i systems necessarily 
involves correlations at different spatial directions this cannot distinguish between classical and quantum correlations 
(not a motivation of 2]), which is the aim of the present paper. 

Here we will show that magnètic susceptibility, when measured along three orthogonal spatial directions, can reveal 
entanglement between individual spins in a sòlid, while magnetization describes their local properties. We also show 
that for a large class of systems magnetization and (zero-field) magnètic susceptibility, when combined in a particular 
way, satisfy a complementarity relation in the quantum-mechanical sense. This macroscopic (thermodynamical) quan- 
tum complementarity relation describes sharing of (quantum) information between entanglement and local properties 
of individual spins in the macroscopic sòlid sample (in an analogy with the relation given above for two qubits). To this 
cnd we will first prové that the sum of magnètic susceptibilities measured along x, y and z directions is a macroscopic 
witness of spin entanglement for a wide class of sòlid state systems. In contrast to infernal energy [Ü EL IS El j the 
present entanglement witness is more general (not only vàlid for special materials 0,0), can be directly measured in 
an experiment, and does not require the complete knowledge of the Hamiltonian of the system. 
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We consider a composite system consisting of N spins of an arbitrary spin length s in a lattice, which is described 
by a spin Hamiltonian Hq. One should mention that entanglement in sòlids may exist in many difFerent degrees 
of freedom, such as spin, spatial or occupation number degrees of freedom. In this letter we will refer to the spin, 
assuming, moreover, that spins can be precisely localized in sites of a lattice. Note, however, that other degrees of 
freedom could also sometimes be represented in the formalism of Pauli spin matrices [l0(. Our results could, therefore, 
also be applicable to these other scenarios. 

In order to study its magnètic response properties, the sòlid is now put in a weak magnètic field, say, directed along 
z-axis and of probe magnitude B p resulting in an additional term H± = B p JZi—i s \ °f the Hamiltonian (here and 
throughout the paper the unit H = 1 is assumed). By s l a {a = x,y,z) we mean here an a-th component of the ï-th spin 
operator in the lattice. Then the Hamiltonian becomes H = Hq + H\. Whcn the system is in its thermal equilibrium 
under a certain temperature T, it is in a thermal state p = e~~ H / kT /Z, where Z = Ti(e~ H ^ kT ) is the partition function, 
and k is the Boltzmann constant. From the partition function one can derive all thermodynamical quantities, e.g. 
the magnetization M z — (l/Z(3)(dZ/dB p ) or the magnètic susceptibility Y?. = (dM z /dB p ), where (3 = l/kT. 

It can be shown that if [H Q , H{\ = 0, magnètic susceptibility is given by 11]: 

l -A 2 {M z ) = ±-{(Ml) - (M z ) 2 ) = — ( V K<; - ( > V ) I . (1) 




kT y ' kT 

where A 2 (M Z ) is variance of the magnetization. 

Microscopically, magnètic susceptibility is, in fact, a sum over all microscopic spin correlation functions {s\s{) for 
the sites i and j. The above is a very important relation as it connects a macroscopic quantity to its microscopic roots 
in the form of the two-site correlation functions. Note, however, that the nonzero value of the correlation function 
does not necessarily imply the existence of entanglement. What we need, loosely speaking, are sufficiently strong 
correlations in all three orthogonal spatial directions and they need to be combined in a specific way to reveal spin 
entanglement. This is the reason why we will now study the sum of magnètic susceptibilities Xx, Xy an d Xz for weak 
probe fields aligned along three orthogonal directions. 

We now show that the expression x.x + Xy + Xz is an entanglement witness. Entanglement witnesses in general 
are observables which (by our convention) have positive expectation vàlues for separable states and negative ones for 
some specific entangled states 01 • The proof is based on the method of entanglement detection using the uncertainty 
relations ^3j- For any separable state of N spins of length l, (for any classical mixture of the products states, each 
appearing with probability w n : p — J2n w nPn ® Pn ® •■■ ® Pn)i one nas 

Ns 

X = Xx+Xy + Xz>j^;- (2) 

We will now prové this inequality for product states of the spins and then the general result will follow 
directly due to convexity of separable states. Note that for an arbitrarily state of spin-s partició one has 
( ( s x) 2 ) +( {s y ) 2 ) +( (s z ) 2 )=s(s + 1) and (s x ) 2 + {s y ) 2 + (s z ) 2 < s 2 . If the thermal state was actually a product 
one of N spins, the variance of magnetization would be the sum of variances of individual spins: x — i^/kT) 
[A 2 (M X ) + A 2 (M y ) + A 2 (M Z )] =(l/kT)Y,i [A 2 (4) + A 2 (sj,) + A 2 (^)] > (l/fcT)(s( s + l)-s 2 ) = Ns/(kT). Note 
that this bound is also vàlid in the general case of separable states due to the convexity of the mixture: 

X = (l/fcT)(A 2 (A4) p + A 2 (M y ) p + A 2 (M z ) p ) 

> (1/fcT) 5>„ Ç [A 2 (4)„ + A 2 ( S ;)„ + A 2 (4)„] 

n i 

> Ns/(kT), (3) 

where index n denotes the n-th subensemble in the mixture, and A 2 {X) p is a variance of an observable X taken in a 
state p. The inequality J2J) is saturated for any pure state which has a maximal magnètic quantum number m with 
respect to any direction (i.e. \j,m= j)). The other extreme case of x — can be achieved for a singlet state of N 
spins where all three variances are equal to zero. 

Therefore, if Xx + Xy + Xz < ^y, the sòlid state system contains entanglement between individual spins. It is 
important to note that all susceptibilities should be taken for zero-fields B p to ensure that they are measured for the 
same thermal state (for the same reason no quantum phase transition at the points of the measurements is assumed) . 
Because the measurement of the magnètic susceptibility has been an experimental routine for long time, it is clear that 
the present approach is an experimentally efficient method for detecting macroscopic spin entanglement. It might be 
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FIG. 1: Detection of entanglement in the xxx Heisenberg spin-i (a) and spin-1 chains (b). The black sòlid curves are the 
theoretical curves from Ref. 0] and represent the temperature dependence of the zero-field magnètic susceptibility x{T) per 
particle in the spin-| (a) and spin-1 chains (b). The red sòlid curves are from our work and represent the right-hand side of 
inequality @. Mathematically, they represent hyperbolae 1/(6T) (a) and 1/(3T) (b). The critical temperatures below which 
entanglement exists in the chains are T c = 1.6J for spins-i and T c = 2J for spins-1 (The data in | l4| were given only for 
T < 2J). 



of particular importance when there is only partial knowledge of systems's Hamiltonian and one thus has to rely more 
on experiment. In what follows we will demonstrate the efncicncy of the method using an exactly solvable model. 

Suppose that the symmetry of the system is such that magnètic susceptibility is equal in all three directions x-x — 
Xy = Xz- This is the case, for example, for the Heisenberg spin lattices with isotropic, but in general inhomogeneuous 
coupling constant Jy-: Hq = . Jijs'íP (here the summation does not need to be constrained to nearest-neighbour 
interactions only). The entanglement criterion now reads as follows: 

X, < (4) 

We apply it to investigate existence of entanglement in the infinite xxx Heisenberg chains of spins-1/2 (Fig. la) 
and spins-1 (Fig. lb), described by a Hamiltonian Hq = J^. SjSÍ+i, at various temperatures. We use the results 
of Ref. [lij where the thermodynamic properties of the Heisenberg spin chains are obtained by the transfer-matrix 
renormalization-group method. The right-hand side of is represented by the red sòlid lines in Fig. ^ The vàlues 
of magnètic susceptibility left to the intersections point of the red and the theoretical curves cannot be explained 
without entanglement. The critical temperatures are T c = 1.6J for spins-^ and T c = 2J for spins-1. 

It is commonly believed that one reconstruets the physics of classical spins in the limit of infinitely large spins 
as suggested by the limit of the commutation relation for normalized spins Si = ^ (i = x,y,z): lim s ^ 00 [S , :c , S y ] = 
lim, s _ >0o -S z — ► 0. Our result shows, however, that the larger the spins are, the higher the critical temperature is 
below which entanglement is present. Perhaps one possible explanation of this effect is that the longer spins produce 
entanglement of higher dimension than the one produced by shorter spins. Thus, for the same coupling strength J 
and temperature the longer spins can develop larger amount of entanglement than the shorter spins, which might 
then persist at higher temperatures. The dependence of an amount of entanglement on dimensionality of subsystems 
was studied in Ref. [ïfij . 

It is interesting to compare the threshold temperature for the spin- ^ chain with the estimates based on other 
macroscopic entanglement witness such as infernal energy. In Ref. [J, |l6( it was shown that the infernal energy can 
reveal the entire bipartite entanglement between neighboring spins as measured by the concurrence . Using these 
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and the results from 01 we obtain that the concurrence vanishes below the threshold temperature of 0.795 J. The 
higher value of 1.6 J as revealed by magnètic susceptibility can be explained by the fact that it, in contrast to internal 
energy, may also detect bipartitivc entanglement between non-neighboring sites and multipartitive entanglemcnt. 

We now turn to the derivation of a macroscopic quantum complementarity relation. We first note that the sum of the 
squares of magnetizations along three orthogonal directions satisfies the relation: (M) 2 = (M x ) 2 + (M y ) 2 + (M z ) 2 < 
N 2 s 2 . This describes the complementarity between properties of individual spins in a sòlid, in analogy with the 
{&x) 2 + (°y) 2 + (cz) 2 < 1 for a single qubit. If one of the observables in the sum, for example (M z ) 2 , takes its 
maximal value of N 2 s 2 (e.g. in state \j — Ns,m = Ns) where j is angular and m magnètic quantum number), the 
other two have to vanish. For the purposes of further discussion we need the following relation between (M) 2 and 



(Ai 2 ) = (M x 2 ) + (M y 2 ) + (M z 2 ): 



(M 2 ) > (^tl) (M) 2 . (5) 

Hcrc follows the proof. Let us denote by | j, m) the joint eigenstates of M 2 with eigenvalues j(j + 1) and M z with 
eigenvalues m. Note that both sides of (JHJ are invariant under rotations in the three-dimensional space. Thus, for any 
given state we can choose such a coordinate system that (M x ) = (M y ) —0 and consequently (M) 2 = (M z ) 2 . Let us now 
dcfinc a new operator K such that K \ j, m) = j \ j, m). Given that pj m are probabilities of finding the system in any 
of the states |j,m), wehave (M z ) = ^2- m pj m m<J2j m Pjmj = (K). To complete the proof we use (M 2 ) — (K{K + 1)) 
and the fact that Ns>j implies Ns(K) > (K 2 ). Thus we have 

m - (^) (Mf >_ (K{K + 1)) - W >_ (^ti) A 2 (A-) > 0. <„ 

We now exploit Eq. and (0 to derive a macroscopic quantum complementarity relation: 

1 kTx (M) 2 



Ns N 2 s 2 



< 1- (7) 



non— local properties local properties 

The left-hand side of inequality Q) can be divided into two parts: Q = 1 — and P = j^] 2 . While P describes 
the local properties of individual spins, Q is associated with quantum correlations between spins in a sòlid. This is 
because Q is proportional to two-site spin correlations for three orthogonal directions (three mutually non-commuting 
observables) and its positive value implies the existence of entanglement (see Eq. J2J). In the extreme case of a 
product state of N spins all aligned along the same direction (e.g. \j — Ns,m = Ns)), their local properties are well 
defined (P = 1) at the expense of no entanglement (Q — 0). In the other extreme case the state of the systems is 
highly entanglcd. Thcn non-local properties are maximally exhibited (Q — 1), at the expense of P = 0. In general, 
the relation (JJJ) describes partial quantum information sharing between local and non-local properties of spins. 

To illustrate the complementarity relation Q we analyze a chain of antiferromagnetically coupled spin pairs - dimers 
- which are themselves uncoupled. This is a correct model for, e.g., Copper Nitrates and many orgànic radicals. The 
Hamiltonian in an external magnètic field of magnitude B is given by 

H = j ° 2j ■ ° 2]+í + B < ■ ( 8 ) 

3 3 

The plot of P, Q and their sum P + Q as a function of the magnètic field B and temperature T is given in Fig. 
For B = 0, the singlet is the ground state and the triplets are the degenerate excited states. For a higher value of 

B, however, the triplet states split and the gap between the singlet and first excited state | ) (a z \ — ) = — | — )) 

decreases. Therefore, in a thermal state at a given temperature as B is increased, the definiteness of non-local 
properties Q decreases because increasingly larger triplet component will be mixed with the sing let HIE!- 0n 
the other hand, as B increases the spins tends to orient themselves all paral·lel to the field, which results in higher 
vàlues of magnctization and thus P, in agreement with the complementarity relation (see Fig. HJl). Increasing T 
generally decreases P + Q as thermal mixing has a destructive character both to Q and P. Note, however, that at all 
temperatures and all vàlues of magnètic field the relation P + Q < 1 is satisfied. 

In conclusions, we show that magnètic susceptibility is an entanglement witness for a wide class of systems. While 
magntetization describes local properties of individual constituents of a sòlid, the magnètic susceptibility specifies 
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FIG. 2: Macroscopic quantum complementarity relation between local (P) and non-local (Q) properties. (a) The plot of 
Q = 1 - 2kTx/N (red), P = 4(M) 2 /7V 2 ( green) and its sum P + Q (black) for a chain of antiferromagnetically coupled spin-i 
pairs (dimers) versus the magnètic field B/J. The temperature is taken to be T = 0.1J. (b) The plot of P + Q as a function of 
magnètic field B/J and temperature T/J. Under all temperatures and vàlues of magnètic field the complementarity relation 
P + Q < 1 is satisfied (see text for further discussion) . 

its spin entanglement. We show that magnetization and magnètic susceptibility satisfy a quantum complementarity 
relation. One of the quantities can thus increase only at the expense of a decrease in the other. This shows quantum 
information sharing in macroscopic quantum systems, such as sòlids. In future, it will be interesting to investigate this 
feature at critical points where (quantum) phase transitions occur. It can be seen from our plot at T = (Fig.Ek) that 
a sudden increase in non-local properties (Q) at the quantum phase transition point B/J = 2 must be accompanying 
with a corresponding sudden decrease in P. Othcrwise, the complementarity relation would be violatcd. 

Our results are not only relevant for fundamental research but also for quantum information science as they give 
the critical vàlues of physical parameters (e.g. the high-temperature limit) above which one cannot harness quantum 
entanglement in condensed matter systems as a resource for quantum information processing. 

C.B. and M.W. were supported by the Austrian Science Foundation (FWF) Project SFB 1506. M.W. is supported 
by the Erwin Schròdinger Institute in Vienna and the Foundation for Polish Science (FNP). C.B. thanks the European 
Commission (RAMBOQ). V. V. thanks European Union and the Engineering and Physical Sciences Research Council 
for financial support. C.B. and V.V. thank the British Council in Àustria. The authors would like to thank B. 
Hiesmayr, A. Fereirra, and J. Kofler on very useful discussions and comments. The work was a part of Austrian- 
Polish Collaboration Programme "Quantum Information and Quantum Communication V" . 



[1] C. Brukner and A. Zeilinger, Phys. Rev. Let. 83, 3354 (1999). 

[2] B. Hiesmayr and V. Vedral, e-print quant-ph/0501015 (2005). 

[3] X. Wang, Phys. Rev. A 66 034302 (2002 ). 

[4] C. Brukner and V. Vedral, e-print quant-ph/0406040 (2004). 



6 



G. Toth, Phys. Rev. A 71, 010301(R) (2005). 

M.R. Dowling, A.C. Doherty, and S.D. Bartlett Phys. Rev. A 70, 062113 (2004). 

L.-A. Wu, S. Bandyopadhyay, M.S. Sarandy, D.A. Lidar, Phys. Rev. A 72, 032309 (2005) . 

S. Ghosh et al, Nature 425, 48 (2003). 

C. Brukner, V. Vedral, and A. Zeilinger, quant-ph/0410138 (2004) (to appear in Phys. Rev A). 
T. Holstein and A.Primakoff, Phys. Rev. Lett. 48, 1098 (1940). 
F. Schwabl, Statistical Mechanïcs, Springer (2002). 

M. Horodecki et al, Phys. Lett. A 223,1 (1996); B. M. Terhal, J. Th. Comp. Sc. 287(1), 313 (2002). M. Lewenstein et al. 
Phys. Rev. A 62, 052310 (2000). 

H. F. Hofmann and S. Takeuchi, Phys. Rev. A 68, 032103 (2003). 

T. Xiang, Phys. Rev. B, 58, 9142 (1998). 

K. Dennison and W. K. Wootters, Phys. Rev. A, 65, 010301, quant-ph/0106058 (2002) and rcfcrcnccs thcrcin. 

X. Wang, P. Zanardi, Phys. Lett. A 301, 1, quant-ph/0202108 (2002), X. Wang, e-print quant-ph/0205049 (2002). 

W. K. Wootters, Quantum Inf. and Comp. 1, 27 (2001). 

M.C. Arnesen, S. Bose and V. Vedral, Phys. Rev. Lett. 87, 01 7901 (2001). 

M. A. Nielsen, PhD Thesis, University of New Mexico (1998), |qüaïü>ph /0011036 



